Towards a semiclassical justification of the 'effective random matrix theory' for 
transport through ballistic chaotic quantum dots 
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The scattering matrix S of a ballistic chaotic cavity is the direct sum of a 'classical' and a 'quan- 
tum' part, which describe the scattering of channels with typical dwell time smaller and larger than 
the Ehrenfest time, respectively. According to the 'effective random matrix theory' of Silvestrov, 
Goorden, and Beenakker [Phys. Rev. Lett. 90, 116801 (2003)], statistical averages involving the 
quantum-mechanical scattering matrix are given by random matrix theory. While this effective ran- 
dom matrix theory is known not to be applicable for quantum interference corrections to transport, 
which appear to subleading order in the number of scattering channels A'^, it is believed to correctly 
describe quantum transport to leading order in N. We here partially verify this belief, by compar- 
ing the predictions of the effective random matrix theory for the ensemble averages of polynomial 
functions of S and 5*^ of degree 2, 4, and 6 to a semiclassical calculation. 

PACS numbers; 73.23.-b,05.45.Mt,05.45.Pq,73.20.Fz 



I. INTRODUCTION 

For electrons moving in a ballistic conductor with 
chaotic classical dynamics there is a minimal time af- 
ter which the wave nature of electrons becomes appar- 
ent: the Ehrenfest timeri The Ehrenfest time rg is the 
time it takes for classical trajectories initially a 'quan- 
tum distance' apart to separate and reach a 'classical 
distance'*^ Here, the quantum distance is the Fermi 
wavelength h/pp, where pp is the Fermi momentum, 
whereas the classical distance can be taken to be the 
system size L or the width W of contacts connecting the 
sample to source and drain reservoirs. For chaotic clas- 
sical dynamics with Lyapunov exponent A, te is then 
determined by the condition W — (?i/pf) exp(ATE), so 
that 



1, PfW 



(1) 



If the Ehrenfest time te is small in comparison to the 
mean dwell time td in the conductor and other time 
scales relevant for quantum transport (such as the de- 
phasing time or the period of an applied AC bias) — 
which is true for most experimentally realized ballistic 
quantum dots^ — , the time threshold it poses for quan- 
tum processes is irrelevant, which explains why quantum 
signatures do not distinguish ballistic conductors from 
their disordered counterparts. Indeed, quantum trans- 
port in ballistic quantum dots with te <SC td and quan- 
tum transport in disordered quantum dots are both de- 
scribed by random matrix theory,^ 

Recently, there has been considerable interest in the 
theoretical question what happens if the Ehrenfest time 
exceeds the mean dwell time td. In this regime signifi- 
cant differences between the manifestations of quantum 
mechanics in ballistic chaotic and disordered conductors 
can occur. It is important to distinguish the effect of 
the Ehrenfest time on quantum phenomena that involve 
the splitting of trajectories only and quantum interfer- 



ence phenomena, which involve the divergence of classi- 
cal trajectories initially a quantum distance apart and 
their joining again. Examples of the former are shot 
j^Qjgg6j7j8j9 g^j^j ^]-^g excitation gap induced by the prox- 
imity to a superconductoriifiiiiiiSiii Examples of the lat- 
ter are weak localizationfi*i^*iSii& universal conductance 
fluctuationSfi2ii2iiS and quantum corrections to the level 
densityi22i2i Here we focus on the first type of phenom- 
ena. The effects of the Ehrenfest time on quantum inter- 
ference effects have proven quite subtle, and will not be 
discussed here (see the references cited above for details) . 

Silvestrov, Goorden, and Beenakker proposed a very 
attractive theoretical picture to describe Ehrenfest-time 
related phenomena in a ballistic chaotic quantum dot.— 
They noted that the dot's classical phase space can be 
divided into a part containing classical trajectories with 
dwell time shorter than te and a part with classical tra- 
jectories with dwell time longer than te. Correspond- 
ingly, the dot's N X N scattering matrix S{e) is written 
as the direct sum of a scattering matrix S'ci(e) for Nc 
'classical channels' and a scattering matrix Sq{e) for A^q 
'quantum channels', with N = Nc + Nn, 
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The Afc-dimensional 'classical scattering matrix' Sc\{£) 
represents fully deterministic scattering from the quan- 
tum dot, where all probability intensity is concentrated 
around one classical trajectory. The scattering phase 
shifts in Sc\{s) follow from the dwell times of the corre- 
sponding classical trajectories. On the other hand, Sq{e) 
represents quantum scattering, where the probability in- 
tensity is divided over a large number of scattering modes 
in all contacts. Writing 



^q(^) 



(3) 



in order to factor out a trivial energy dependence arising 
from the fact that all trajectories contribution to Sq have 
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a minimal dwell time TEr^*^ Silvestrov, Goorden, and 
Beenakker proposed that the statistical distribution of 
Sq{s) is that of random matrix theoryfi^ provided A'q is 
large. 

Since the original proposal of Ref. 12, it has been mi- 
derstood that this 'effective random matrix theory' does 
not provide a faithful description of all signatures of 
quantum transport. For example, the effective random 
matrix theory predicts that the weak localization cor- 
rection to the conductance of a ballistic quantum dot is 
Ehrenfest-time independent, whereas both microscopic 
semiclassical theory and numerical simulations find an 
exponential dependence on te/td Weak local- 
ization is a quantum interference effect, which arises as a 
correction to subleading order in the total channel num- 
ber TV. The effective random matrix theory has been suc- 
cessful in predicting and explaining Ehrenfest-time de- 
pendences that appear to leading order in N, such as shot 
noise, ^■^■'^■^ the proximity-induced gap in a quantum dot 
coupled to a superconductor ^^"i^^i^^i^'^i^^ the density of 
transmission eigenvalues f24 and the probability distribu- 
tion of lifetimes of quasibound states in chaotic quantum 
dots.^ 

Whereas the effective random matrix theory has been 
compared to the results of accurate numerical simula- 
tions in all four examples mentioned above, the micro- 
scopic verification of the effective random matrix theory 
is limited to a theoretical construction of the decomposi- 
tion |j2J and to the comparison of the effective random 
matrix theory and microscopic calculations of the te de- 
pendence of the shot noise power of a chaotic quantum 
dot-^ and the density of states of a quantum dot cou- 
pled to a superconductor (an 'Andreev quantum dot') for 
energies much larger than the proximity-induced energy 
gappii (For energies comparable to the energy gap, the 
theory of Ref. [ij makes use of an ansatz that is similar 
in spirit to the ansatz of effective random matrix theory.) 
The comparison to a microscopic calculation of shot noise 
amounts to a test of the effective random matrix theory 
for traces of polynomial functions of the scattering ma- 
trix S and its hermitian conjugate of degree 4. The 
comparison to the density of states of an Andreev quan- 
tum dot at high energies addresses of the statistics of 
polynomials of Sci{s)S*i{—e) of arbitrary degree, but not 
of Sq. The aim of this article is to provide the next step 
towards a microscopic verification of the effective ran- 
dom matrix theory by comparing a semiclassical theory 
of the ensemble average of a trace of a degree-six poly- 
nomial function of S and with the predictions of the 
effective random matrix theory. Since the hypothesis of 
the effective random matrix was formulated after the mi- 
croscopic calculations of the Ehrenfest-time dependence 
of shot noise and the density of states in an Andreev 
quantum dot, we believe that this calculation is the first 
nontrivial test of the effective random matrix theory. 

In Sec. |n] below we review the predictions of the ef- 
fective random matrix theory for the averages of traces 
of polynomials of S and of degree two, four, and six. 
We follow with a semiclassical calculation of these aver- 
ages in Sec. IIIII and conclude in Sec. IIVI Details of the 
calculations can be found in the two appendices. 



II. PREDICTIONS OF THE EFFECTIVE 
RANDOM MATRIX THEORY 

In our calculation, we consider a ballistic quantum dot 
coupled to electron reservoirs through ballistic point con- 
tacts. Hence, the scattering matrix 5* acquires a block 
structure S = Sij, where the indices i and j label the 
point contacts. The dimension of the block Sij is Ni x Nj , 
where Ni is the number of channels in the ith point con- 
tact. Each block Sij can be decomposed into a 'classical' 
and a 'quantum' scattering matrix as in Eq. (j^J. 

We are interested in averages of the form 



Q2 



^(tr5,,(£i)5|,(£2)), 



N 

^{trS^^{s^)Sye2)S,iie,)SUe,)), 
^{trS,^{e^)Sl{s2)S^iie3) . . . Sl{ee)). (4) 



The polynomials Q2 and Q4 describe, e.g., the con- 
ductance and shot noise power of a quantum dot with 
normal-metal contactSf2& whereas polynomials of higher 
degree are necessary for a theory of transport and equi- 
librium properties of quantum dots with superconduct- 
ing contactsii^ According to the effective random matrix 
theory, these averages have the structure 



Or: 



-_F„(l,...,n)TE/TD 



,-F„(l,...,n)rE/rDgRMT 



where the function Fn is defined as 



i^„(l,...,n) = l + ^5](-l)^ 



(5) 



(6) 



and Q'^' and Q^^""" are obtained from Eq. Q by replacing 
S by Sc\ and Sq and N by Nc and A^q, respectively. 
The classical averages QJj' readi^ 



Ar2F2(l,2)- 

N,Nj^^kSJl 



Ar2F4(l,2,3,4)' 

cl _ NjNjSjkSjiSimSjn 

^ ~ Ar2F6(l,2,3,4,5,6)- 



(7) 



The random matrix average of the trace of a degree two 
polynomial of S and 5^ is equal to the classical average^Si 



QT^ = Qf- 



(8) 



However, the higher-order random matrix averages are 
different. For Q4 one has2& 
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V4 



iV3F2(l, 2)^^2(3, 4) iV3F2(3, 2)^2(1, 4) 
whereas the random matrix average Qq reada^ 



N,N,NkNiF4{l,2,3A) 
NiF2{l, 2)F2{3, 2)F2(3, 4)^^2(1, 4) ' 



(9) 



^RMT 

We 



N,N,NkNiNrnNnFe{l,2, 3, 4, 5, 6) 
iV4F2(l, 2)^2(3,4)^^2(5, 6) ' 7V4i^2(3,2)i^2(5,4)F2(l,6) ^2(1, 2)^2(3, 2)i^2(3, 4)^2(5, 4)^^2(5, 6)^2(1, 6) 
N,N,NiN^S,kSin , 7V,Ar,7VfcA^,A^™iV„F4(l, 2, 3, 4)^4(1,4,5,6) 



+ 



N^F2(1, 2)F2(3, 6)F2(5, 4) ^2(1, 2)^2(3, 2)^2(3, 4)^2(5, 4)^2(5, 6)^2(1, 6)^2(1, 4) 
N,N,NiNraNnFii3, 4, 5, 6)d,k N,NjNkN„,N„Fi{l, 4, 5, 6),5,/ 



7V5^^2(1, 2)^^2(3, 4)^2(5,4)^^2(5,6)^2(3, 6) 7V5^^2(3, 2)^^2(1, 4)^2(5,4)^^2(5,6)^2(1, 6) 



+ .. 



(10) 



where the dots . . . refer to the simultaneous cychc per- 
mutations (l,2,i,j) — > (3,4, fc,/) — > (5,6, m,n). 

Together, Eqs. |(SJ, 0, ©, ©, and specify the 
prediction of the effective random matrix theory for the 
averages of the traces of polynomials of the scattering 
matrix and its hermitian conjugate for polynomials up 
to degree 6. 



lively, and m and n label the propagating modes in these 
leads. The widths of the entrance and exit contacts are 
Wj and Wi, respectively. The trajectory a connects the 
entrance contact j to the exit contact i. The components 
p± and p'^ of its momentum perpendicular to the lead 
axis upon entrance and exit, respectively, are compatible 
with that of the modes n and m in the corresponding 
leads, 



III. SEMICLASSICAL CALCULATION 

We now perform a semiclassical calculation of the av- 
erages Q2, Q4, and Qe defined in Eq. Q above and show 
that they agree with the predictions of the effective ran- 
dom matrix theory as described in the previous section. 
Although semiclassical c alcul ations of Q2 and Q4 exist, 
see Refs. ■27.30..31i3^ and l6l9ll3M respectively, we briefly 
review their derivation in order to establish the context 
for the calculation of Qg. (We are not aware of a semi- 
classical calculation of the energy dependence of Q4, how- 
ever.) 

In our calculation, we take the limit h while keep- 
ing the ratios te/td and Ni/N , and the products Toei/h 
fixed. The latter condition implies that the functions Fn 
defined in Eq. © remain constant in the limiting pro- 
cedure. Note that the channel numbers Ni diverge in 
this limit. The divergence of the channel numbers does 
not affect Q2, Qi, or Qq, however, because these depend 
on the ratios Ni/N only. Since te oc ln(l/?i), see Eq. 

the condition that the ratio te/td is kept constant 
in the limiting procedure implies that the dwell time td 
diverges as well. The divergence of td removes any de- 
pendence on the non-universal short-time dynamics in 
the quantum dot. 

Starting point of our calculation is an expression of the 
dot's scattering matrix S' as a sum over classical trajec- 
tories a )?9i?ii?^ 



{Sij)r 



f ■kTi 



1/2 



(11) 



where i and j label the exit and entrance leads, respec- 



p± = ±whn/Wj, n = 1, . 
p'^ = -inrhm/Wi, m — 1, 



(12) 



(Here and in the remainder of this article, primed vari- 
ables refer to the exit contact.) Further, Sa is the clas- 
sical action of trajectory a and is its stability ampli- 
tude. The latter is defined as 



dy 



'1/2 



(13) 



where y is the coordinate perpendicular to the axis of the 
entrance contact, see Fig. ^ and the partial derivative 
is taken at constant p±. The classical action Sa is a 
fimction of p± and p'^. The spatial coordinates y and 
y' of the trajectory a upon entrance and exit can be 
expressed as derivatives of Sa , 



y - 



dSa 

dp± 



y 



dSa 

'dp'^ 



(14) 



For simplicity of notation, the Maslov index and other 
phase shifts are included in Sa ■ 

We are interested in the scattering matrix Sij at dif- 
ferent values of the energy e. The energy e enters in Eq. 
()ll|l through the energy dependence of the action Sa , 



Sa{£) = Saisa) + (e - £0)^0 



(15) 



where Eq is a reference energy and ta is the duration of 
the trajectory a. We neglect the energy dependence of 
the stability amplitudes Aa- 

Using Eq. ifTT)) . the traces Q„, n = 2,4:, 6, are expressed 
as double, quadruple, and sixfold summations over clas- 
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Pi. cavity 



FIG. 1: Schematic drawing of one of the contacts to the 
quantum dot. The coordinate y as weU as the perpendicu- 
lar component of the momentum p± refers to the direction 
perpendicular to the lead axis. In the sum 1)11^ over classi- 
cal trajectories there is no restriction on y, but only discrete 
values of p± are considered. 



sical trajectories, 



12 = 



2W,Wj 



J2 AiA2e'^^'- 



S2)/h 



(16) 



with similar expressions for Qi and Qe- Here and in 
the remainder of this article we use the indices 1 and 2 
to denote ai and a2 if possible without confusion. The 
transverse momentum components p± and p'^ at the en- 
trance and exit contacts of the n trajectories involved in 
the trajectory sum for Qn satisfy 



P±,2 



P'±.3^ 



P±,n-1 = P±.n, 
I I 

P± 



(17) 



where n — 2,4, 64S In the summation (|16|l . the magni- 
tudes of all transverse momentum components are taken 
equal to the quantized values l(T^ . However, for large 
channel numbers we may replace the summation over the 
quantized momenta by integrations over p_L and p_L' , so 
that 



Q2 - ^ / dp±,i'^P'-L,i E MA2e^^^--^^^"\ (18) 

ai,Q2 



again with similar expressions for Q4 and Qg- In Eq- 
and its equivalents for (54 and Qq, one integrates over 
the perpendicular momentum components p±_ and of 
the trajectories ai, . . . ,q;„_i with odd indices. The ini- 
tial and final transverse momentum components of the 
classical trajectories q;2, . . . , q;„ with even indices are de- 
termined by the conditions (|17|l . 

For the calculations below, we assume that the classi- 
cal dynamics inside the quantum dot is uniformly hy- 
perbolic. This assumption simplifies the calculations, 
although it is believed not to affect the final results in 
the universal limit td ^ L/vp^^'^^^ With uniformly 
hyperbolic classical dynamics, the optimal phase space 



coordinates for a Poincare surface of section taken in 
the interior of the quantum dot are coordinates s and 
u taken along the stable and unstable directions in phase 
space. In view of this, we replace the trajectory sum (|18|l . 
which is taken over trajectories with specified transverse 
momentum components p± and p'^ at entrance and exit 
contacts, by a trajectory sum over trajectories with spec- 
ified stable and unstable phase space coordinates s and 
u' at entrance and exit contacts, respectively. Referring 
to ApD.IXIfor details, we find that the traces Q2, Q4, and 
Qq can also be expressed as 



-S2)/n 



(19) 



{2T:hy 

X J2 Ai...Aie'^^^---^^'^/'', (20) 
/ dsidu\dsy,duLdst^du'c: 

{2TThY J i i o 5 



i-...-Se)/h 



(21) 



where the classical trajectories a^, fi = 1, . . . , n satisfy 
the conditions 



■ 1 -^n— 1 — 3n , 



(22) 



with n = 2, 4, 6. Further, in Eqs. the stability 

amplitudes are defined as 



Aa — 



du' 



du 



-1/2 



(23) 



and the classical actions Sa{s,u') are Legendre trans- 
forms of the original classical actions Sa(pA^,p'j^), so that 



ds 



du' 



(24) 



For the interpretation of Eqs. ((T^ " (|^ . one should 
keep in mind that the phase space coordinates s and u 
are defined only locally. That means that the coordinate 
transformation (a;,p_L) — > (s,m) can only be made for 
pairs of trajectories that enter or exit the quantum dot 
at nearby phase space points. This poses no problems for 
our calculation, because only classical trajectories that 
exit the quantum dot at nearby positions and with close 
momenta have an action difference AiS that varies suffi- 
ciently slowly as a function of the phase space coordinates 
to give a finite contribution to Qn, n — 2,4,6. Indeed, 
the explicit calculations of Q2, Q4, and Qe in the follow- 
ing subsections show that the entire contribution to Qn 
comes from trajectories for which the differences of the 
stable or unstable phase space coordinates are small. 
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A. Calculation of Q2 

The leading contribution to Q2 arises from equal tra- 
jectories ai = a2- In that case, the action difference 



AS = Si- S2 
= (£1 - £2)*, 



(25) 



where t is the common duration of the trajectories ai and 
a2- The stability amplitudes are equal, A\ = A2 — A. 
The factor = in Eq. (fT^ provides the Jaco- 

bian necessary to replace the integration over the unsta- 
ble phase space coordinate u' at the exit contact by an 
integration over the unstable phase space coordinate u 
at the entrance contact i22*^ii2^ This way, Q2 is expressed 
in terms of an integration over the phase space coordi- 
nates s and u at a Poincare surface of section taken in 
the entrance contact. Equivalently, the double integral 
over s and u' can be represented by a phase space inte- 
gral over any Poincare surface of section taken anywhere 
along the classical trajectory a\ — a2^ Following the 
latter strategy, we write Q2 as 



dq / dtidt: 



P,ih)Pj{t2)e'^''-'''^''''+''^/'' 
' 2nh{ti +t2)N ' 

(26) 

where q refers to the phase space coordinate at which the 
reference surface of section is taken, Pi{ti) selects only 
those q for which the classical propagation of q ends up 
in contact i after time ti, and Pj{t2) selects q for which 
the classical propagation of the time-reversed of q ends 
at contact j after a time t2- We divided by t = ii -f i2 
to cancel a spurious contribution from the freedom to 
choose the reference surface of section at an arbitrary 
point along the trajectory ai = a2- Replacing Pi and Pj 
by classical probabilities to reach the contacts i and j for 
an arbitrary phase space point g. 



-*/^°, P,{t) = 



N 



(27) 



the integration over q contributes the total phase space 
volume 2'k?iNty) of the quantum dot and we find the 
well-known result 



Q2 = 



(28) 



N^F2{1,2) 
upon integration over ti and t2- 

B. Calculation of Q4 

The trace is expressed as a quadruple sum over 
classical trajectories, see Eq. (jSOJ. The typical config- 
uration of four trajectories that contributes to Q4 was 
pointed out in Refs. I^IsS The trajectories ai, 02, as, 
and a4, which are paired close to entrance and exit, have 
a small-angle encounter at which the pairing of the tra- 
jectories is interchanged, see Fig. |5| Before arriving at 




t ,^ t 




FIG. 2: Left: Four classical trajectories that give the random 
matrix contribution to the average Q4. These trajectories 
have a small-angle encounter which fully resides inside the 
quantum dot. Right: Schematic drawing of the encounter 
together with the definitions of the various times used in the 
text. The encounter region is shown thick. 



the encounter, ai and a2, and and are paired. Af- 
ter the encounter, a2 and 013, and 014 and ai are paired. 
The encounter may reside fully inside the quantum dot, 
as in Fig. 121 so that the pairs of trajectories are uncorre- 
lated when they exit and enter the quantum dot, or the 
small-angle encounter may touch one or two of the lead 
openings, as in Figs. |21 and 2) respectively, so that exit 
or entrance of the four trajectories is correlated. While 
the importance of such small-angle encounters of classi- 
cal trajectories was first realized for weak localizationjiiSi 
they play a crucial role in all semiclassical theories of 
quantum transport 

Since all four trajectories are fixed once we have spec- 
ified ai and 03, we need to sum over ai and 0:3 only. 
In order to parameterize ai and we take a Poincare 
surface of section chosen at an arbitrary point during the 
encounter. We then use the coordinate g of a reference 
point at the Poincare surface of section, as well as the 
differences S3 — si and U3 — wi of the stable and unsta- 
ble phase space coordinates to parameterize ai and asi^ 
The trajectories ai and a2, and 03 and a4 originate from 
the same contact and with the same stable phase stable 
space coordinate, hence 



Sl 



S2, S3 = S4. 



Similarly, 



U4 = Ui, U2 = U3. 



(29) 



(30) 



Since the trajectory pairs have the same stable and un- 
stable phase space coordinates at the entrance and exit 
contacts, respectively, the conditions and hold 
for encounters that do not touch the contacts as well as 
for encounters that touch the contacts. 

We first discuss the case of Fig. in which the en- 
counter resides fully inside the quantum dot. We closely 
follow Ref. 33, in which the same configuration of tra- 
jectories was considered in the limit te/td — > 0. The 
encounter region is defined as the segment of the trajec- 
tories for which |si — S3I < c and \ui — U3I < c, where 
c {ppWy^'^ is a classical cut-off below which the clas- 
sical dynamics can be linearized. The precise choice of 
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W is irrelevant in the classical limit, see the discussion 
following Eq. (|37|l . Hence, the duration of the encounter 



^cnc \ I / 

A \[U3 



ui)(s3 - Sl)| ' 



(31) 



where A is the Lyapunov exponent for the classical dy- 
namics in the quantum dot. We parameterize the dura- 
tions , , , and ta^ of the four trajectories involved 
using ionc and the durations ti2, ^23, ^34, and ^41 of the 
four stretches connecting the encounter region with the 
lead openings. 



tcnc + tl2 + til, 
tcnc + ti2 + ^23, 
^cnc + ^34 + ^23, 
^cnc +^34 +^41- 



(32) 



Schematically, the definitions of tone and of the times ii2, 
^23, ^34, and til are shown in the right panel of Fig. |21 
The action difference AS is the symplectic area enclosed 
by the four trajectories,'^^'^^ plus a contribution from the 
energy differences, 

AS = (S3 - Si)(u3 - ui) 

+ tcnc(£l — S2 + £3 ^ £4) 

+ tuiei - £2) + t23(e3 - £2) 

+ i34(e3 - £4) + t4i(ei - £4)- (33) 

Note that the enclosed phase space areas are conserved 
along the motion of the trajectories. In particular, this 
means that the action difference AS is independent of 
where the Poincare surface of section is chosen. 

Integrating over the position of the Poincare surface of 
section, we then find 



dh2dt23dt3iUiPj{h2)Pkit23)Plit34)Bdtil) 



X / d{u3 - ui)d{s3 - si) 



Toe 



;A5/;i-tc„c/TD 



2Trfitr 



(34) 



In this equation, the factor tone in the denominator can- 
cels a spurious contribution to the integral from the free- 
dom to choose the Poincare surface of section anywhere 
along the encounter region. The classical probabilities 
Pi{t4i), Pj{ti2), Pk{t23), and Pi{t34) are defined as in 
Eq. l|?7|l . The factor exp(— tone/To) is the probabihty 
that the trajectories do not exit the quantum dot during 
the encounter stretch. (If they do, the encounter touches 
the lead opening. This case is treated separately below.) 

Taking U3 — ui to be positive (while adding a factor 2 
in Eq. (jHU), we perform the variable change 



U3-ui = c/a, S3 - si 



(35) 



With the new integration variables, the integration do- 
main is — 1 < X < 1 and 1 < a < l/\x\. Further, 




FIG. 3: Schematic drawing of four classical trajectories that 
give the random matrix contribution with one Kronecker delta 
to the average Q4. These trajectories have a small-angle en- 
counter touches one of the lead openings. The Kronecker 
delta appears for the entrance lead index in (a) and for the 
exit lead index in (b). 



tcnc = A^^ ln(l/|a::|). Integrating over ti2, t23, t34, and 
t4i, and (7, and performing a partial integration to x, we 
then find 



N,N,NkNiF4il,2,3,A) 
' mF2{l,2)F2{3,2)F2{3A)F2{l,4:) 
< /^xx^4i,2,3.4)A^2sinM^ (3g) 

Jo 



where we omitted a term that is an oscillating function of 
r = /Tim. the limit ?i — > 0. The remaining integral over 
X can be evaluated in the limit ?i ^ at fixed te/td, 
and we arrive at the final result 



Q 







N,N,NkNiFi{l,2,i,A) 
>4i;^2(l, 2)^2(3,2)^2(8,4)^^2(1, 4) 



X e 



-teF4(1,2,3,4)/td 



(37) 



where the Ehrenfest time te is defined as, cf. Eq. 



TE 



A 



In?' 



1 



(38) 



It is important to notice that the precise value of the 
phase-space cut-off c enters only through the Ehrenfest 
time Te- However, since the Ehrenfest time appears in 
the combination te/td only and since td — > 00 in the 
classical limit taken here, any change of c by a factor of 
order unity will not affect the final result. It is because of 
this that we did not need to carefully specify the classical 
length scale W entering into the cut-off c. 

If the encounter region touches one of the lead open- 
ings, as in Figs. |2t and b, a slight modification of the 
above calculation is called for. For dcfiniteness, we con- 
sider an encounter that touches the entrance lead open- 
ing, as shown in Fig. This implies j = I, so that 
the corresponding contribution to Q4 is proportional to 
5ji. The action difference AS is given by Eq. (|33|l with 
ti2 = ^34 = 0. Since the trajectories ai and a3 must be 
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correlated upon entry for the configuration shown in Fig. 
the encounter time tcnc is bounded by 



1 



^ c 1 c 
< ionc - Y In : < - hi : 

A |u3-ui| A |s3-si| 



(39) 



Then, proceeding as in the calculation of and inte- 
grating over ^23 and tn , we find 



'4 N^F2{3,2)F2{1A) 

d{si - S3)d{ui - U3) 



4 

dt^YicPj (^cnc) 



X 



2'Khtp 



y. Qi{s3—Sl)(u3—Ul)/h + itoac{ei— 62+63— ei))/h (^4Q^ 

In order to perform the integrations in Eq. (|40|l . we take 
U3 — ui to be positive and perform the variable change 



U3 — Ul = — , — Si — CX(T. 

(T 



(41) 



With the new integration variables, the integration do- 
main is— 1 < X < 1, 1 < a < 6'''*™% and < icnc < 
A^^ ln(l/|a;|). Hence, with r = /h and integrating over 
tcnc and a, the integral ipn|l becomes 



7V3F2(3,2)F2(1,4) 

^ ^^ 2sin(rx) ^j._^(;^ 2,3,4)/ArD 



-teF4(1,2,3,4)/td 



Ar3F2(3,2)F2(l,4)' 



, (42) 



where, again, we omitted terms proportional to sinr. 
Similarly, for one finds 



Q 



H) 



iV3F2(3,4)F2(l,2) 



-TE-F4(l,2,3,4)/7-D 



(43) 



Finally, if the encounter touches both lead openings, 
one has 



Q 



gl 



TD / rfti(ii2 7 difcOj! 



X / d(s3 - Si)d(w3 - ui)- 



2nh ' 



(44) 



where and Pj are classical probabilities given in Eq. 
(|77|l . and ^2 are the propagation times between the 
surface of section and the exit contact and the entrance 
contact and the surface of section, respectively, and 
^enc = ^1-1-^2- The action difference AS is given by 
Eq. with ti2 = t23 = ^34 = ^41 = 0. In order to 
ensure that ai and as are correlated upon entrance as 
well as exit, we require 



1 1 c 
<i < - In : 

A |S3-5l| 



1 



i2 < T In 



A |u3-ui|' 



(45) 




FIG. 4: Schematic drawing of four classical trajectories with a 
small-angle encounter that touches both lead openings. Such 
configurations of classical trajectories give the classical con- 
tribution to Q4. 



Performing integrations to ti — ^2 and 1*3 — ui and ab- 
breviating X = (s3 — si)c, one finds 



-At 

dx 2 sin(a;r) 



63-6i) 



(46) 



The integral over x has to be performed in the limit r 
00 at fixed ratio te/^. Since 



dx 2 sin(a;r) 



X TT 



9{te - t) 



(47) 



if r ^ 00, where 9{x) — 1 if x > and otherwise, we 
find 



NiNjS.kSjl _ g_rEF4(l,2,3,4)/rD)^ (48) 



iV2P4(l,2,3,4)" 



Together, Eqs. (jSH), 1321), 63), and (gHIl reproduce the 
prediction of the effective random matrix theory. 

We note that there is a close connection between the 
appearance of Kronecker deltas involving the contact in- 
dices i, j, fc, and I, and the configurations of classical 
trajectories that contribute to Q^: each encounter region 
that touches a lead opening gives rise to a Kronecker 
delta and, conversely, each Kronecker delta derived from 
an encounter that touched a lead opening. Hence, the 
four contributions to Q4 — the third line in Eq. |(7J) and 
the three terms in Eq. — , each of which have a dif- 
ferent product of Kronecker deltas, are in a one-to-one 
correspondence with the four configurations of classical 
trajectories shown in Figs. |21 Et, Et*, and^ 

We should point out that, without dependence on the 
energy arguments, i.e., after setting £1 = £2 = £3 = £4, 
the semiclassical calculation of Q4 was first performed 
by Whitney and Jacquod.^ The same results can also be 
inferred from an earlier calculation by Agam, Aleiner, 
and Larkin, who used a different formalism.^ The energy 
dependence of the four terms we calculate here is non- 
trivial, however, and clearly reveals the structure of the 
classical trajectories underlying the four different terms 
in the average. 
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FIG. 5: Schematic drawing of six classical trajectories that 
contribute to Qe without encounters that touch the lead open- 
ings. The configuration shown in (a) consists of two sep- 
arate two-encounters. There are two more configurations 
of this type, which are obtained by the cyclic permutations 
(f,2) — > (3,4) (5,6). The configuration shown in (b) con- 
sists of one three-encounter. 

El 

While our calculation of closely followed that of 
Braun et al.^ our calculation of the remaining three con- 
tributions to Q4 differs considerably from that of Ref . ^3§. 
Braun et al. do not consider encounters that touch the 
lead openings. Instead, they calculate the remaining con- 
tributions to Q4 using the 'diagonal approximation' for 
the trajectory sums. Although this approximation gives 
the correct result if te ^ td, it is based on a different 
class of trajectories than the ones we considered (see also 
Ref. 0). In the diagonal approximation, one considers 
the case that the four trajectories ai, 02 , as, and are 
pairwise equal (rather than pairwise close), but without 
imposing any relation between the two pairs. For exam- 
ple, one admits trajectories ai = a2 and — 0:4 where 
a2 and exit with the same transverse momentum p± 
but at a classically different spatial coordinate y. When 
summed over the full family of trajectories, such configu- 
rations appear with rapidly oscillating phases, and their 
net contribution vanishes These rapidly oscillating 
phases disappear only if all four trajectories pass through 
the contact at equal positions and angles (up to quantum 
uncertainties), as is the case for the configurations shown 
in Figs. 12^1 and b. 



C. Calculation of Qe 

The generic configuration of classical trajectories that 
contributes to Qq is shown in Fig. [S] Instead of dealing 



with trajectories with one small-angle encounter, one now 
has to consider configurations of trajectories with two 
small-angle encounters or of trajectories with a 'three- 
encounter A three-encounter arises when two en- 
counters overlap, so that all six trajectories are within a 
phase space distance c simultaneously. In order to avoid 
confusion, we refer to non-overlapping encounters of two 
pairs of trajectories as 'two-encounters' for the remain- 
der of this section. We first consider the case that all 
encounters fully reside inside the quantum dot. The case 
that encounters touch the lead openings will be discussed 
afterwards. 



Without encounters that touch the lead openings, one 
distinguishes two types of contributions to Qq: the con- 
tribution from the case in which the trajectories undergo 
two separate two-encounters, and the contribution from 
the case in which the six classical trajectories involved 
undergo a three-encounter. These are shown in Fig. 
and|3l3, respectively. The configuration shown in Fig.[S^ 
shows only one configuration with two two-encounters. 
There are two more configurations, which are obtained 
by the cyclic permutations (I,2,«,j) (3,4, fc, Z) — > 
(5, 6, TO, n). 



Because the two two-encounters in Fig. do not 

overlap, their contribution Qg to Qe factorizes. 
The stretch connecting the trajectories ai and 02 
to their joint entrance contact j contributes a fac- 
tor Nj /NF2{1,2); the other stretches connecting the 
encounters to the contacts contribute similar factors 



to 



Further, using the results of the pre- 



vious subsection, we find that the two-encounter of 
the trajectories ai, q;2, aa, and 0:4 in Fig. con- 
tributes a factor -^4(1, 2, 3, 4) exp(-TE^4(l, 2, 3, 4)/rD). 
Similarly, the two-encounter involving the trajec- 
tories Q!i, 0:4, 0:5, and ae contributes a factor 
-F4(l,4,5,6)exp(-rE^4(l,4,5,6)/rD). Finally, the 
stretch of the trajectories ai and that connects 
the two two-encounters in Fig. contributes a factor 
l/i^2(l,4). Using that i^4(l,2,3,4) i^4(i4,5,6) ^ 

^6(1,2,3,4,5,6) + ^2(1, 4), we find that is given 

by 



iV,iVjiVfc7V/7V™7V„e-^EF6(l,2,3,4,5,6)/rD 

7V6F2(1, 2)^^2(3, 2)^^2(3, 4)F2(5, 4)^2(5, 6)^2(1, 6) 



f4(l,2,3,4)f4(l,4,5,6) _.,;^,(i,4)/.o 
i^2(l,4) 



, (49) 



where the dots . . . represent two more terms obtained by the cyclic permutations (1, 2, i, j) (3, 4, k, I) —>■ (5, 6, to, n). 
As before, the superscript jSJi) refers to the figure that shows the corresponding configuration of classical trajectories. 

Although the calculation of is similar in spirit, it proves rather cumbersome. Referring to App.^for details. 
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(e) (f) (g) 

FIG. 6: Schematic drawing of six classical trajectories that provide the remaining contributions to Qe- The configuration 
shown in (a) consists of two separate two-encounters, where one two-encounter touches a lead opening. The configuration 
shown in (b) consists of one three-encounter, which partially touches one lead opening. The configuration of (c) consists of a 
three-encounter that fully touches one of the lead openings. The configuration of (d) has two two-encounters that each touch 
lead openings. The configuration of (e) has a three-encounter that partially touches two lead openings. The configuration of 
(f) has a three-encounter that partially touches one, and fully touches the other lead opening. Finally, the configuration of (g) 
fully touches both lead openings. The contributions from (a) and (b) have one Kronecker delta involving the lead indices, the 
contributions from (c), (d), and (e) have two Kronecker deltas, the contribution of (f) is zero, and the contribution of (g) has 
four Kronecker deltas. 



we find 



'6 7V6^2(1, 2)F2(3, 2)F2{3, 4)F2(5, 4)^2(5, 6)^2(1, 6) 

" 7V6F2(1, 2)^2(3,2)^2(3,4)^^2(5, 4)^2(5,6)^2(1, 6) 



f4(l,2,3,4)j^4(l,4,5,6) 
^^2(1, 4) 

^6(1,2,3,4,5,6). 



I _ g-rEF2(l,4)/rD 



(50) 



Combining these two contributions, we find that the terms proportional to exp(— 2rE/rD) cancel, so that the final 
result is uniformly proportional to exp(— te/td). 



|5j ^ jV,jVj-iVfciV;iV„,jV„e-"^^«(l-^'3,4,5,6)/ro 

6 " iV6F2(l,2)F2(3,4)F2(5,6)i^2(l,6)F2(3,2)F2(5,4) 

^4(1,2,3,4)^4(1,6,5,4) ^ F4(3, 4, 5, 6)^4(3,2,1,6) , ^^4(5, 6, 1, 2)i;^4(5,4, 3, 2) 



^6(1,2,3,4,5,6) 



(51) 



In order to compute the contribution to Qq of encounters that touch the lead openings, we have to consider the 
configurations of classical trajectories shown in Fig. El FigureEt show trajectories with two two-encounters where one 
encounter touches the lead opening. Figure IHb shows a configuration of classical trajectories with a three-encounter 
where the three-encounter partially touches a lead opening. In both cases, the escape of only one pair of trajectories is 
correlated, so that the contribution of the configurations of Figs.lHt and b has one Kronecker delta. In addition to the 
configurations shown in Fig. and b, there five more contributions to Qq that are obtained by cyclic permutations 
of the trajectory labels and by the exchange of entrance and exit contacts. As before, the contribution from the 
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configuration with two two-encounters factorizes, 

A^,iVziV„e-"^^«(l'2.3.4,5,6)/ro 



iV3F2(l,2)F2(3,4)F2(5,6) 

Af,Affc7V„e-"E^«(l'2,3,4,5,6)/ro 

7V3^^2(3, 2)^2(5,4)^2(1, 6) 



jV,iVfcJ4(l,2,3,4)(5„ 
iV2i?^2(2, 3)^^2(1, 4) 

N,NiFi{3A,5,2)S,n 
iV2F2(3, 4)^^2(5, 2) 



g-TEF2(l,4)/TD _^ 

teF2{5,2)/to I 



(52) 



Details of the calculation of the three-encounter are again left to the appendix. The result is 



7V,iViiV„e-"^^«(l'2,3.4,5,6)/rD 
' iV3f^2(l, 2)^2(3,4)^^2(5, 6) 



N,NkFiil,2,3,A)Sun A 
7V2^^2(2,3)F2(1,4) V 

iVj-jV;f4(3,4,5,2)^j-„ / 
7V2f^2(3, 4)^^2(5, 2) V 



-te-F2(1,4)/td 



-TE-F2(5,2)/rD 



7V3F2(3, 2)^2(5,4)^2(1, 6) 

Combining both contributions, we again find that the terms proportional to exp(— 2te/td) cancel, so that 



(53) 



N,NjNkNiNn5^rnFi{l,2, 3, 4) 



+ 



N,NjNkNiN„,6jnF43A,5,^) 



N^F2{l,2)F2{3,2)F2{3,4)F2i5,6)F2ilA) N^F2{3,2)F2i3,A)F2i5,4)F2{l,6)F2i5,2) 



+ . . 



X e 



^rEF6(l,2,3,4,5,6)/TD 



(54) 



There are two different types of contributions to Qq with two Kronecker deltas: Two Kronecker deltas that both 
involve lead indices of entrance (or exit) leads, or a product of two Kronecker deltas, where one involves lead indices 
of the entrance lead and one involves lead indices of the exit lead. Only a configuration of trajectories with a three- 
encounter contribute to the former, see Fig. EJ;, whereas configurations of trajectories with a three-encounter as well 
as configurations of trajectories with two two-encounters contribute to the latter, see Fig. EJl and e. Referring to the 
appendix for calculational details regarding the configurations of trajectories with a three-encounter, here we simply 
quote the results, 



Q, 



Etc) 



N,N,NiNnS,kS 



N,N,NkNraSilS, 



7V4F2(1, 2)^^2(3,4)^2(5, 6) 
N.N^NiN^S^kSin 



7V4i^2(3,2)F2(5,4)i^2(l,6) 
N,NJNkNJ,„^SJl 



-teFs (1,2, 3,4,5, 6)/rD 



N.NjNkNidjnSkr. 



N^F2{1, 2)F2{3, 6)F2(5, 4) N^F2{1, 4)^2(3, 2)^^2(5, 6) N^F2{1, 6)F2{3, 4)^2(5, 2) 



X e 



-TEFe(l,2,3,4,5,6)/TD 



(55) 



(56) 



Finally, we have to consider the case that there is cor- 
related escape of one pair and one triplet of trajectories 
and correlated escape of two triplets of trajectories, see 
Figs.Ef and g, respectively. One verifies that there is no 
contribution to Qe for the configuration shown in Fig. |Sf . 
For the configuration shown in Fig. |Ht one finds 

nEl) _ N^NjSikSjmSjlSjn , -rEF6(l,2, 3,4,5. 6)/tdN 

iV2F6(l,2,3,4,5,6)^' ^• 

(57) 

Together, Eqs. (jSIl, (Ell, lESJ, ISEl), and combine 
to the effective random matrix theory ansatz for Qq. 



IV. CONCLUSION 

In this article we calculated the Ehrenfest-time depen- 
dence of the ensemble averages of polynomial functions 
Qn of the Nx N scattering matrix S* of a ballistic chaotic 



quantum dot and its hermitian conjugate , for polyno- 
mials of degree n = 2, 4, and 6. For all cases our cal- 
culations of the leading-in-A^ averages agreed with the 
effective random matrix theory of Silvestrov, Goorden, 
and Beenakker.^2 

The detailed semiclassical calculations of Q2 and Q4 
were included in this article because they are a prerequi- 
site for the semiclassical calculations of Qq. For the cal- 
culation of Q4 a few guiding principles would have been 
sufficient, however. These are: (i) the classical trajec- 
tories contributing to Qi have at most one small-angle 
encounter, (ii) each encounter lasts one Ehrenfest time 
te, and (iii) Q4 becomes equal to the classical limit if 
te is much larger than the mean dwell time td , whereas 
Q4 equals the random matrix limit Q^^"^ if te <IC td. 
From the first two guiding principles, both of which are 
well established in the literaturefi*^ and the exponential 
distribution of classical dwell times, one concludes that 
Qi is of the form 

g4 = A + Bexp(-TE/rD). (58) 
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The third guiding principle then fully determines the co- 
efficients A and B and, hence, Q4. 

No such shortcut exists for the calculation of Qq. The 
relevant configurations of semiclassical trajectories have 
up to two encounters, hence we expect a general form 



Qe = A + B cxp{-te/to) + C cxp{-2te/tb). (59) 



The detailed calculation of this article was needed to 
show that the coefficient C = 0. In this sense, the semi- 
classical calculation of Qq provides the first nontrivial 
test of the effective random matrix theory. 

The fact that the coefficient C in Eq. H59(l is zero is at 
the heart of the effective random matrix theory ansatz: 
It implies that (i) all trajectories with dwell times longer 
than te are treated equally and (ii) that te is the only 
time marking the boundary between classical and quan- 
tum effects. Whereas our calculation indirectly verifies 
these claims for leading-in-A^ averages of low-degree poly- 
nomials of S and 5^, they are not true for subleading- 
in-A^ averages. For subleading-in-A^ averages, more than 
one integer multiple of te may determine quantum ef- 
fects. Known examples are weak localization, to which 
only trajectories with a dwell time larger than 2te con- 
tribute, although the weak localization correction is sup- 
pressed cx exp(— te/td) with an exponent containing te, 
not 2TEr*i^*^ the leading quantum correction to the 
spectral form factor K{u>) of a quantum dot with broken 
time-reversal symmetry^Si which has oscillations propor- 
tional to exp(2ia;TE), exp(3icijrE), and exp(4ia;TE), con- 
ductance fluctuations in a bulk chaotic conductor (rather 
than a chaotic quantum dot),^^ and the mean square cur- 
rent pumped through a chaotic quantum dot with time- 
dependent potentials. 

Having verified the effective random matrix theory 
ansatz for leading-in-A^ averages for polynomials in S and 
up to degree six, one wonders whether this verifica- 
tion can be extended to polynomials of arbitrary degree. 
In our calculation, the simple result predicted by the ef- 
fective random matrix theory was obtained only after a 
tedious calculation showing that the coefficient C in Eq. 
(|59|l above is zero. The possibility of the appearance of 
a term oc exp(— 2te/td) caused a large increase in com- 
plexity upon going from polynomials of degree four to 
polynomials of degree six. We have searched for a sim- 
plification in our calculation that would allow us to cancel 
all contributions oc exp(— 2rE/TD) at an earlier point in 
the calculation. However, since the terms proportional to 
exp(— 2te/td) were obtained from qualitatively different 
configurations of classical trajectories, we have not been 
able to find such a simplification. Without such a sim- 
plification, the task of extending the present calculation 
to polynomials of arbitrary degree seems too daunting to 
accomplish. 




FIG. 7: Poincare surface of section of lead opening, together 
with the standard coordinates y, p±, as well as the phase space 
coordinates s, u, taken along the stable and unstable direc- 
tions in phase space. The trajectory pair ai, 02, which have 
equal perpendicular momentum components px,i = p±,2 are 
mapped to a trajectory pair ai, 02 with equal stable phase 
space coordinates, §1 = S2- The unstable phase space coordi- 
nates do not change in the mapping. 
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APPENDIX A: TRAJECTORY SUM WITH 
FIXED COORDINATES s AND u' 

In this appendix we show that one can replace the stan- 
dard formulation in which the traces Qn are expressed in 
terms of integrals over classical trajectories with specified 
transverse momentum components at entrance and exit 
contacts by an integral over classical trajectories with 
specified stable and unstable phase space coordinates at 
entrance and exit contacts, respectively. 

We first consider a single integral over the transverse 
momentum p± in the entrance contact, 

(AI) 

where the trajectory has transverse momentum pj_ at 
its entrance contact, 

P±.i=P±.2=P±, (A2) 

and transverse momentum ^ at its exit contact, /i = 
1,2. The stability amplitude is given by Eq. (|13|l . The 
classical action is a function of p± and p^, fi — 1, 2. 

In order to rewrite K in terms of an integral over clas- 
sical trajectories with a fixed stable phase space coordi- 
nate s upon entrance into the quantum dot, we consider 
a Poincare surface of section taken at the entrance con- 
tact, see Fig. [71 There are two sets of coordinates to label 
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trajectories piercing this Poincare surface of section: the 
spatial and momentum coordinates y and p^, and the 
stable and unstable phase space coordinates s and u, see 
Fig. [7| The coordinate axes for s and u do not need to be 
perpendicular, but the coordinates are normalized such 
that dsdu — dydp±^ The phase space points where ai 
and a2 pierce the Poincare surface of section are indi- 
cated in the figure. 

To each pair of two sufficiently close classical trajec- 
tories ai, a2 we now assign another pair ai, where 
the trajectory ci^ is obtained from by slightly chang- 
ing the boundary condition in the entrance contact while 
keeping the boundary condition upon exit the same. This 
implies that the unstable phase space coordinates of 



and 



are equal, fi — 1,2. In order to uniquely de- 



fine the pair ai, a2, we require that ai and 0:2 have the 
same stable phase space coordinate s upon entry, and 
that P± I — P±.i — P±,2 — P± 2- construction of ai 
and a2 is shown in Fig. [71 

For this construction, one verifies that the action dif- 
ference of the original and primed trajectory pair is not 
changed, 



5*1 — S2 — Si — 5*2, 



(A3) 



provided we replace the action S, which is a function of 
p±, by the action S, which is a function of the stable 
phase space coordinate s, obtained by Legendre trans- 
form of the original action S. Also, using the normaliza- 
tion dsdu = dp±dy and Eq. (|A2|) . one finds 



a three-encounter that fully resides inside the quantum 
dot. 




FIG. 8: Definition of the four measures of the encounter 
duration used in the text: tcnc is the total duration of the 
three-encounter, t^nc is the duration of the part of the en- 
counter that all trajectories involved in the three-encounter 
are within a phase-space distance c, and ts and are the 
durations of the encounter stretches where one pair of trajec- 
tories has already diverged from the remaining two pairs. 



In order to calculate , we take a Poincare surface 
of section at a point that all six trajectories are within a 
phase space distance c. We parameterize the trajectories 
using phase space coordinates s and u along the stable 
and unstable directions in phase space. At the Poincare 
surface of section, one has 



Sl = S2, S3 = S4, S5 = Sq 
U2 = U3, Ui = Us, U6 = ui. 



(Bl) 



dpi 



dp'j_ 



dyi 



^1/2 



dy2 



-1/2 



where 



1 ds 




-1/2 


9p'l,2 


-1/2 


dui 










s = 


SJ = 


S~2- 





(A4) 



(A5) 



The same procedure can be applied to integrations over 
the transverse momentum p'^ in the exit contacts. In this 
case, one replaces the integral over by an integral over 
the unstable phase space coordinate u' , and replaces pairs 
of trajectories with equal p'^ by pairs of trajectories with 
equal u' . The result of this program is precisely Eqs. 
(^J-ljll]) of the main text. 



(These equalities continue to hold if the three-encounter 
touches the lead opening.) The action difference AS 



AS 



(U3 - Mi)(s3 - 
-|- Eiti — £2^2 



Sl) + ("5 ~ Ul)(s5 - S3) 

^...-sete- (B2) 



In order to parameterize the durations of the six tra- 
jectories involved, we separate each duration tf^, p, = 
1, 2, . . . , 6, into segments before and after the encounter 
and one segment of duration i^^cnc inside the encounter 
region. Similarly, we define 



APPENDIX B: THREE-ENCOUNTER OF 
CLASSICAL TRAJECTORIES 

1. Three-encounter fully inside the quantum dot 



A5o„c = ("3 - Wl)(s3 - Sl) -f (U5 - Ui)(s5 - S3) 

+ £1^1, one — £2^2,0110 + . . . — £6^6,0110- (B3) 



In this appendix we describe the details of the cal- 



culation of Q 



This is the contribution to from 



Integrating over the time segments outside the encounter, 
we then arrive at 
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N,N,NkNiN„,Nn 



7V6F2(l,2)F2(3,4)F2(5,6)i^2(l,6)F2(3,2)F2(5,4) 



with 



d{s3 - Sl)d{s5 ~ Sl)diu3 - Ul)d{u5 - Ml) jASo„c/?»-tc„c/rD 



Here tcnc is the total duration of the encounter (the time that at least two trajectories are close together), 



tc. 



A niin(|s3 - si|, |s5 - si|, |s3 - S5I) A niin(|M3 - |m5 - |w3 - U5I) ' 
and is the time that all three trajectories are close together, see Fig. |H1 

c 



cnc \ /I II I 

A max(|s3 - si|, |S5 - si|, |S3 



A max(|M3 - |u5 - |u3 - U5I) 



(B4) 



(B5) 



(B6) 



(B7) 



The factor t^j^^ in the denominator cancels a spurious contribution to the integral arising from the freedom to choose 
the Poincare surface of section at any point during three-encounter. The integration domain in Eq. ljB5p is given by 
the conditions max(|it3 — |m5 — |m3 — M5I) < c and max(|s3 — si|, \s^ — si|, |s3 — S5I) < c. 

We rewrite the integral (|B5|I using s = max(|s3 — si|, |s5 — si|, |s3 — S5I), s' — min(|s3 — si|, |s5 — si|, |s3 — S5I, 
u = max(|M3 — ui\, |it5 — Mil, |u3 — U5I) and u' = min(|u3 — mi|, |u5 — ui|, |u3 — U5I) as our integration variables. After 
some tedious algebra, one arrives at 



TD 



dsdu 



s/2 



ds' 



u/2 



du'- 



-t^„^-F6(l,2,3,4,5,6)/rD 



g-F4(l,2,5,6)t,/rD-F4(l,6,3,2)t„/rD ^ g-F4(l,2,3,4)t,/rD-F4(l,6,3,2)t„/rD _|_ 



us' + u's — US 



us + u's' — u's 



h h 

-F4(l,2,3,4)t,/TD--F4(5,4,l,6)t„/rD , 



us + u's' — us' 

COS h COS ■ 



cos ■ 



COS ■ 



(B8) 



where the dots indicate cyclic permutations (1,2) 
(3,4) -> (5,6), and with 



/ 

cnc ~ 


1 , c 
-In - 4 
A u 


- - In-, 
A s 




ts - 


A 1^7' 

1 u 
A^V- 






tu = 




(B9) 



The total duration of the encounter is t, 



cnc — "-one 



We then perform the variable change 



s ~ c/(T, s' ~ cx' /a, 
u — cy<7, u' — cyy'a. 



(BIO) 



With these new variables, the integration over a can 
be done and cancels the factor t'^^^ in the denominator. 
Writing r = c^/fi., the remaining integral then reads 



f'^' dx'dy' t ydyy^«(i>2'3.4,5,6)/Aro 

Jo Jo 

X I "(■2;')f'4(l,2,5,6)/ATD(-y/-)F4(l,6,3,2)/ArD __|_ (^^/^F4(1,2,3,4)/Atd (^y/^F4(l,6,3.2)/ArD __|_ ^ 

X [cos{yr{x' + y' — 1)) + cos(?/r(l + x'y' — y')) + cos(?/r(l + x'y' — x')) + cos(yr(l — x'y'))] 
(^/)F4(i,2,3,4)/Aro(y/)F4(5,4,i,6)/Aro ^ j [cos(?;r(a;' + y' - x' y') + cos(yr(x' - ?;'))]}■ 



(Bll) 
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In order to evaluate Eq. (|B11|) . we write each product {x'Y^{y'Y^ as 

{x'YHyT = (1/2)^^+^^ + [{xT - {i/2r]{i/2r- + [{yT - (1/2)^^(1/2)^^ + [{x'r - (1/2)'^] [(2/')^^ - (1/2)^1 

(B12) 

and evaluate each of the four terms in Eq. ljB12|) separately. In Eq. (|B12|) . the exponents ei and 62 represent the 
exponents ^4(1, 2, 5, 6)/Atd etc. in Eq. IjBlip . The separation in Eq. (|B12|) makes sense, because both ei and £2 are 
sent to zero if the classical limit ?i ^ at fixed ratio te/td is taken. Hence, the first term in Eq. ljB12|l approaches 
unity, whereas the other terms are nonzero only if x' or y' are small, or both. 

The four terms in Eq. l|B12p are evaluated separately. For the first term in Eq. (|B12I) , the integrals to x' and y' can 
be done exactly. Omitting a prefactor {l/2y^~^'^^ , which is sent to 1 in the classical limit taken here, and abbreviating 
ry = Fg{l, 2, 3, 4, 5, 6) / Ar^, we find that the integrals from the first term give 

h = ^^-^ — / dx'dy' [ dyyy''[cos{yr{x' + y' - I)) + cos{yr{l + x'y' - y')) 



+ cos(j/r(l + x'y' — x')) + cos(j/r(l x'y')) + cos{yr{x' + y' — x'y')) + cos{yr{x' — y'))] 
^'^^^ ^ dyy''-^ h9c{ry) cos(ry) + g^iry) sin(ry)] , 

oy 



(B13) 



where 



/ ^ I J 1-cosw r sinw /m.N 
gdv) = / dw , g^iv) = / dw . (B14) 



Then, performing a partial integration to y and omitting terms that are oscillating with r, we find 

h = ^^!2!Z f [g^ry) cos{ry) - gsiry) sm{ry)] , (B15) 

Ja y 

In this expression we can take the limits r ^ 00 and 77 0, keeping r^^ — e^'^E^6(i,2,3,4,5,6)/rD jjxed. We then find 

h = -F6(l,2,3,4,5,6)e-"^^«(i'2'3.4,5,6)/rD_ (B^g) 

For the second term in Eq. ljB12|) . we first perform the integration to y' , with the result 

' dx'dy' dyyy'^[{xy -{l/2Y][cos{yr{x' +y' -l))+cos{yr{l+x'y' ^y')) 



+ cos{yr{l + x'y' — x')) + cos(yr(l — x'y')) + cos{yr{x' + y' — x'y')) + cos{yr{x' — y'))] + . . . 
= ^ J^'^\,' J\yy.[(^,'y _ (^^,).j| ^^cos(,r.O -cos(ry) +jl^-.Ocos(yKl -^0) ^ _ ^^^^^ 

where we abbreviated e = ^4(1, 2, 3, 4)/ArD. Performing a partial integration to y, omitting oscillating terms, we then 
find 



»l/2 



, dx'ix'^-am- , 



^^6(1,2,3,4,5,6) 



'{x" - (1/2)') r ^^y" ^' cosjryx') - cosjry) + (1 - a;')cos(yr(l - x')] 

Jo y x'{i-x') 

(BI8) 



The integral in the first line is of order 1/r and can be neglected in the limit r 00. The integral in the second line 
is of order e cx I/Xtu, so that it can also be neglected. This is best seen by evaluating the integral after expanding 
{x')"^ — (1/2)*^ « eln(2a::') and replacing y'' by 1, which gives a finite value of the integral. Hence, we conclude that, in 
the classical limit r ^ 00 at fixed r^'^/^'^o ^ have 

I2 = 0. (B19) 
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Similarly, one finds that the third term in Eq. ljB12|) does not contribute to /. 

For the remaining contribution I4, we note that the only nonzero contribution can come from small x' and small 
y'. Neglecting x' and y' with respect to unity, one finds that the contribution from the first term between brackets 
{. . .} in Eq. IjBlip gives zero. In the same approximation, after two partial integrations the second term gives 



4ATDeie2 f^^^ dx'dy' 



x'y' Jo 



dyy^ {x'Y^ {y'Y^ sm{yrx') sm{yry') + 



(B20) 



where ei = ^4(1, 2, 3, 4)/Atd and €2 ~ ^4(5,4,1,6) (or cyclic permutations), and 77 — i^6(l, 2, 3, 4, 5, 6)/Atd. We then 
shift variables x' = x" /y and y' — y" /y and integrate over y. 



h = 



4ATDeie2 



dx"dy" 



x"y" 



{x'T{y'T[{^B.^{x"y)) 



7r2(ei +62-77) 

f4(l, 2, 3, 4)F2(5, 4, 1, 6) .,f,(i,2.3,4.5,61/.n 



{1/2^-"^-"^] sin(rx") sh\{ry") + 



-TEF2(l,4)/rD 



(B21) 



where we used that ^4(1, 2, 3, 4) + ^4(5,4, 1,6) - 
^6(1,2,3,4,5,6) = ^2(1, 4). As before, the dots . ..refer 
to the two additional terms obtained by cyclic permu- 
tation (1,2) — > (3,4) (5,6) of the expression written 
above. 



The final result for Qg 
h +h into Eq. (|B4)| . 



is found by substituting / = 



2. Three-encounter that touches the lead openings 

If the three-encounter touches the lead openings, es- 
cape of the trajectories involved in the encounter is 
no longer uncorrelated. In order to deal with three- 
encounters that touch the lead openings, the initial part 
of the calculation of the previous subsection has to be 
modified. However, the final integrals will all be of the 
form of the integrals h, I2, and I4 calculated above. Here 
we discuss the configurations of Figs.Eb and c in detail. 
The calculation of the remaining configurations (Figs.EJ?- 
g) proceeds in a similar manner. In the discussion here 
we further limit ourselves to the labeling of trajectories as 
shown in the figure. For Fig. Eh, this means that we con- 
sider the exit of trajectories ai and to be correlated, 
whereas the exit of trajectory as, as well as all entrances 
are uncorrelated. There are five more contributions to 
Qe- two of these arise from configurations where the exit 
of Qfi and as or the exit of 03 and a^ is correlated; three 
more arise from configurations with correlated entry of 
two trajectories, and uncorrelated exits. For Fig. Et, we 
take the exits of all trajectories to be correlated, whereas 
the entrances are not. There is another contribution ob- 
tained by reversing the roles of entrance and exit. 

In order to describe a three-encounter that touches 
the lead opening as shown in Fig. ^jp, we again take a 
Poincare surface of section at a point that all six trajec- 
tories are within a phase space distance c. In the configu- 
ration of Fig.Eb, this 'central' part of the three-encounter 
does not touch a lead opening. We ensure that as has 
moved away from the other trajectories before arriving 




FIG. 9: Definition of the four measures of the encounter 
duration used in the text: tcnc is the total duration of the 
three-encounter, f^nc is the duration of the part of the en- 
counter that all trajectories involved in the three-encounter 
are within a phase-space distance c, and ts and t' are the 
durations of the encounter stretches where one pair of tra- 
jectories has already diverged from the remaining two pairs. 
The right end of the encounter is at the lead opening. 



at the lead opening by requiring 

min(|ui - U3I, |u5 - U3I) > |ui - U5I. (B22) 
The escape of ai and is correlated if 

0<t' < tu, (B23) 

where 



1 max(lu3 - Mil, |m3 - u^l) 

tu = -r In ■ ■ , 

A U5 - wi 



(B24) 



and t' is the duration of the stretch of the encounter 
immediately adjacent to the lead opening, during which 
Qfi and 015 are correlated with each other, but not with 
03 , see Fig. El We then find 



M) 



il 



iV5F2(l,2)F2(3, 4)^^2(5, 6)F2(3, 2) 



/, (B25) 



where 



I = I d{s3 - si)d{s5 - si)d{u3 - ui)d{u5 ~ ui) 

dt' e'^-Scno-tonc/rD 



(B26) 
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with 



The action difference AiSo„c is given by Eq. IjBSp above. 



^^'^ " A^'Vax(|u3 - 

+ t1^ — ^ 

A max(|si 
1 , max(|,si — 
A mm Si 



(B27) 



c 



- Sal, - |S3 - S5I 



Sal: pi ^ Ssli P3 ^ S5 



(B28) 



Rewriting the integral in terms of the variables s = 
max(|si -S3I, |si -S5I, |s3 -S5I), s' = min(|si -53], |si - 
S5I, |s3 - S5I), u = max(|M3 - ui\, \u3 - U5I), and u' = 
|u5 — and integrating over t' , we arrive at 



;/2 



u/2 



-tonc-P'6(l,2,3,4,5,6)/rD r 



1 _ g-*..-F4(l,4,5,6)/rD 



X COS ■ 



■ COS ■ 



+ 



2g-tsFi{l,2.3,4)/rD 



COS ■ 



us' + u's — us' 



We now change variables 



s = c/a, s' — cx' /a, u = cya, u' = cyy'cr. 



(B29) 



(B30) 



Writing r = /Ti and integrating over a, the integral then reads 



f'\y' r ^j^j^i.e(1.2.3.4.5.6)A^[l_(^,)i..(1.4.5.6)A. 

7r^F4(5,4, 1,6) Jo Jq Jo 



X |[(2;')-P4(1,2,5,6)/Atd _|_ (-^/^F4(3,4,5,6)/ArDj 

X [cos(yr(l + y' x' — y')) + cos(y7'(l — x' — y') + cos{yr{l — x'y')) + cos{yr{\ — x' + y'x'))] 

+ 2(a;')^"^''''^'*^/^"°[cos(yr(a;' + y' - x'y')) + cos{yr{x' - y'))]} • (B31) 

The remainder of the calculation is identical to that of the previous subsection: We write 

(x'r = l + ((a;')^-(l/2r), (B32) 



where e represents ^4(1, 2, 5, 6)/ArD, ^4(8, 4, 5, 6)/ArD, or i^4(l, 2, 3, 4)/ArD, and evaluate the two terms in Eq. (|B32|I 
separately. The calculation of the first term is identical to that of I2, see Eqs. (jB17t - (jB19|l . and vanishes. The 
calculation of the second term is identical to that of I4, see Eq. (|B21|I above. We thus find 



^ -^4(1, 2, 3, 4) ^^;;Fe(1.2,3.4,5.6)/rD 
i^2(l,4) 



I _ g--rE-F2(l,4)/TD 



(B33) 



r 



Substitution of Eq. (|B33I) into Eq. (|B25|) and addition of where 
the remaining five configurations of trajectories obtained 
by relabeling or by interchanging entrance and exit gives 
Eq. @. 

Proceeding in a similar manner for the case shown in 
Fig- Et I we find j _ 



N,N,NiNJ,k5,, 



7V4F2(l,2)F2(3,4)i^2(5,6) ' 



(B34) 



TD / d{s3 - si)d{s5 - si)d{u3 - Ui)d{u5 - Ui) 



dt'^ 



piASo„c/?i-tcnc/rD 



(B35) 
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with 

t - t' 

We require 



1 inax(|si - S3I, |5i - S5U53 - S5I) 
A min(|si - S3I, \si - S5I, \s3 - S5I) ' 

(B36) 



, 1 c 

^ ^cnc ^ ~\~ 71 i i i i TT 

A max(|si - S3I, |si - S5I, |S3 - S5I) 
+ - In 



A max(|ui - U3I, |ui - usl, |u3 - U5I) ' 

(B37) 



in order to ensure that the three-encounter indeed 
touches the lead opening at a point that aU trajectories 
involved in the encounter are separated by phase space 
distances less than the cut-ofF c. The integration over t'^^^ 
can be done immediately. Rewriting the integral in terms 
of the variables s — max(|si — S3I, |si — S5I, |s3 — S5I), s' — 
min(|si -S3I, |si -S5I, |s3 -S5I), w = max{|Mi -M3I, \ui - 
M5I, IM3-U5I), and u' = min(|ui-U3|, {ui-u^l, IM3-M5I), 
and performing the variable change 



c/cr, 



cx' /a, 



cya, u' = cyy'a, (B38) 



we arrive at the integral 



2ATDr2 



1/2 



dx'dy' / ydy[l - yMlM5,6)/Xrn^ r(^/)F4(l,2,5,6)Aro ^ 



^2^^6(1,2,3,4,5,6; JO JO 

X [cos(yr{l + y x — y')) + cos{yr{l — x — y') + cos(?;r(l — x'y')) + cos(yr(l — x' + y'x')) 
+ cos{yr{x + y' — x'y')) + cos(yr(a;' — y'))]- 



(B39) 



In this integral, the contribution from the term 1 in Substituting this result into Eq. (|B34p and adding what 

[1 — y^s(^'2>3>4,5,6)/ATDj vanishes. Writing {x'Y — (1/2)'^+ one obtains after interchanging entrance and exit leads, 

{{x'Y — (1/2)'^), with e = ^4(1, 2, 5, 6)/ArD etc., we arrive one arrives at Eq. ifS^ . 
at integrals identical to the integrals Ii and I2 considered 

in the previous subsection. We thus find The remaining configurations of trajectories with a 

three-encounter (Figs.Et; andlHJl are treated in the same 

/ = g-rEi=^6(i,2,3,4,5,6)/rD^ (B4q^ manner. 
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